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The  results  yield  reasonably  accurate  short-term  rain  rates  for  nearly  98 
percent  of  the  rain  period.  Either  model  readily  determines  the  percent  of 
time  the  average  rain  rate  or  any  selected  rain  rate  is  equalled  or  exceeded. 

Of  significance  is  the  fact  thpt  the  mean  daily,  hourly,  or  instantaneous 
rain  intensity  for  any  duration  tends  to  be  greatly  exceeded  for  at  least 
5  percent  of  the  time.  The  graphs  and  computer  programs  given  not  only 
facilitate  the  rapid  estimation  \of  short-term  rainfall  for  almost  any 
situation,  but  also  serve  to  improve  understanding  of  short-term  rainfall 
spectra. 


PREFACE 


This  report  on  the  determination  of  short-term  rainfall  spectra  for  a  broad  range  of  rainfall  regimes 
is  part  of  an  Army- wide  investigation  of  Battlefield  Obscuration.  The  results  are  directed  toward 
the  design  engineer,  who  is  concerned  with  the  impact  of  rainfall  on  materiel,  and  also  toward  the 
field  commander,  whose  ground  operations  may  be  hindered  or  thwarted  because  of  rain.  Intense 
rain  causes  soil  erosion  or  flooding,  hampers  mobility,  penetrates  the  protective  covering  of  electrical 
equipment,  reduces  visibility,  causes  microwave  attenuation,  and  contributes  to  die  structural  weakening 
or  collapse  of  roofs.  The  rainfall  models  presented,  either  general  or  explicit  for  a  given  average 
rain  rate,  yield  quick  estimates  of  daily,  hourly,  or  instantaneous  rainfall  for  much  of  the  respective 
rainy  period.  The  required  parameters  are  the  total  rainfall  and  the  duration  of  rain  in  the  appropriate 
units  of  time.  Total  rain  and  the  number  of  days  of  rain  per  month  or  year  are  usually  available  in 
most  climatic  summaries.  An  hourly  rainfall  index  developed  during  this  study  affords  a  means  of 
approximating  the  duration  of  rain  in  hours  per  year  if  not  known.  The  derivation  of  instantaneous 
rainfall  is  limited  to  single  storms.  The  graphs,  equations,  and  computer  programs  included  in  this 
report  facilitate  the  determination  of  the  frequency  of  a  selected  rain  rate  or,  given  a  selected  frequency, 
the  rain  rate  that  would  be  equalled  or  exceeded  per  given  interval  of  time. 

The  work  was  accomplished  under  project  QG652C0D01,  Work  Unit  01,  “Frequency  and 
Distribution  of  Natural  Battlefield  Obscurants.  “Appreciation  is  extended  to  all  personnel  in  the 
Battlefield  Environmental  Effects  Group  and  to  Bruce  Zimmerman  in  the  Intelligent  Systems  Group 
of  the  Geographic  Sciences  Laboratory,  U.S.  Army  Engineer  Topographic  Laboratories,  for  their 
helpful  suggestions. 

The  work  was  performed  under  the  supervision  of  Dr.  Donald  W.  Dery,  Branch  Chief,  Battlefield 
Environmental  Effects  Group;  Regis  J.  Orsinger,  Chief,  Land  Combat  Systems  Division;  and  Bruce 
K.  Opitz,  Director,  Geographic  Sciences  Laboratory.  COL  Alan  L  Laubscher,  CE,  was  Commander 
and  Director,  and  Walter  E.  Boge  was  Technical  Director  of  the  U.S.  Army  Engineer  Topographic 
Laboratories  during  the  report  preparation. 
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SIMPLE  ANALYTICAL  METHODS  FOR  ESTIMATING  SHORT-TERM  RAINFALL 


INTRODUCTION 


Information  on  short-term  rainfall  interests  agriculturists,  hydrologists,  geomorphologists,  and 
construction  or  radar  engineers.  Excessive  rain  can  cause  flooding,  decrease  of  soil  traction,  erosion 
of  terrain,  structural  weakening  or  collapse  of  roofs,  and  malfunction  of  jet  engines  because  of  water 
ingestion.  Intense  rain  may  penetrate  the  protective  covering  of  electrical  equipment,  reduce  visibility, 
and  cause  microwave  attenuation.  Military  planners  are  concerned  with  the  impact  of  rain  on  mobility 
and  trafficability.  Data  for  rain  accumulation  over  short  periods  of  time  are  generally  not  available. 

Routine  precipitation  data  for  much  of  the  world  consist  of  the  average  monthly  and  annual  total 
amounts  and  the  corresponding  number  of  days  with  measurable  rain  (an  amount  which  varies  from 
one  country  to  another,  as  greater  than  or  equal  to  .25  mm  or  .  1  mm  for  the  U.S.  or  Southeast  Asia, 
respectively).  However,  a  number  of  weather  service  bulletins  list  daily,  hourly,  or  shorter-term 
precipitation  rates,  particularly  for  excessive  rain,  for  various  durations,  including  5-,  10-,  and 
15-minute  periods  (U.S.  National  Oceanic  and  Atmospheric  Administration  1983;  Panama  Canal 
Company  1971). 1 2  Except  for  instantaneous  data,  which  require  more  sensitive  rain  gauges,  all  hourly 
rainfall  and  most  daily  rainfall  utilized  refer  to  rain  accumulations  equal  to  or  greater  than  .25  millimeter. 
For  those  stations  in  the  tropics  for  which  the  0.1  mm  limit  applies,  the  number  of  days  of  rain  is 
not  very  different  from  that  based  on  the  .25  millimeter  amount  (the  difference  is  of  significance 
in  the  high  latitudes,  where  the  rainfall  is  much  gentler). 

From  selected  weather  bulletins  or  else  from  specialized  experiments  at  a  variety  of  sites,  scientists 
have  developed  a  number  of  instantaneous  precipitation  models  (Bussey  1950;  Briggs  and  Harker 
1969;  Davis  and  McMorrowl976;  Crane  1980;  Tattelman  and  Scharr  1983;  Jones  and  Wedland  1984; 
Bertheland  Plank  1985). 3  4  5  67  8  Recently,  members  of  the  U.S.  Air  Force  Geophysical  Laboratory 
(AFGL)  reviewed  in  detail  and  evaluated  several  models  (Tattelman  and  Grantham  1982).9  Certain 
of  these  models  derive  short-term  rainfall  (hourly  or  instantaneous)  directly  from  total  rain  and  the 
number  of  days  of  rain  per  year  (month).  In  some  instances,  additional  parameters  such  as  temperature 
and/or  a  moisture  index  are  also  employed.  The  AFGL  model  (Tattelman  and  Scharr  1983)  consists 
of  a  set  of  linear  regression  equations,  each  specific  for  a  given  percent  frequency  per  year  that  an 
estimated  instantaneous  rainfall  is  equalled  or  exceeded. 


’ti  S  National  Oceanic  and  Atmospheric  Administration  (NOAA),  Climatological  Data,  annual  summary  (since  1950). 

2 Panama  Canal  Company,  Climatological  Data,  Canal  Zone  and  Panama,  annual  reports  1957-1971  [Balboa  Heights,  Canal  Zone). 
-’Howard  t.  Bussey,  ", Microwave  Attenuation  Statistics  Estimated  from  Rainfall  and  Water  Vapor  Statistics, "  Proceedings  of  the  Institute 
of  Radio  Engineers,  vol.  38  (July  1950):  781-785. 

4Allen  R  Davis  and  Daniel  J.  McMorrow,  Stochastic  Models  for  Deriving  Instantaneous  Precipitation  Rate  Distributions,  Air  Weather 
Service  (Scott  AFB,  Illinois:  U.S.  Air  Force  Environmental  Technical  Applications  Center,  1976,  AWS-TR-76-263). 

Crane,  " Prediction  of  Attenuation  by  Rain,"  IEEE  Transactions  on  Communications.  COM-28,  9  (1980):  1717-1733. 
“Douglas  M.A,  Jones  and  Wayne  M.  Wedland,  "Some  Statistics  of  Instantaneous  Precipitation. "  Journal  of  Climate  and  Applied  Meteorology 
(1984):  1273-1285. 

7J.  Briggs  and  J  A.  Harker,  '  "Estimation  of  the  Duration  of  Short-Period  Rainfall  Rates  Based  on  Clock-Hour  Values. '  The  Meteorological 
Magazine  97  (London,  1969):  246-252. 

8Robcr,  0.  Berthel  and  Vemon  G.  Plank,  Time  Durations  of  Rain  Rates  Exceeding  Specified  Thresholds  (Hanscom  AFB,  Massachusetts: 
U.S  A,r  Force  Geophysics  Laboratory,  1985,  AFGL-TR-0122). 

9Paul  Tattelman  and  D.D.  Grantham.  "A  Survey  of  Techniques  for  Determining  Short-Duration  Precipitation  Rate  Statistics  "Air  Force 
Surveys  in  Geophysics,  no.  441  (Hanscom  AFB,  Massachusetts:  U.S.  Air  Force  Geophysics  Laboratory.  1982,  AFGL-TR  82-0357). 


Berthel  and  Plank  (1985) 10  showed  that  percent  frequencies  of  instantaneous  rainfall  could  be 
estimated  by  the  employment  of  the  coefficient  of  variation  of  the  rain  rate.  This  coefficient  was 
found  to  be  a  function  of  the  ratio  of  maximum  instantaneous  rainfall  divided  by  mean  instantaneous 
rainfall.  Their  technique  requires  that  both  parameters  for  the  ratio  be  obtained  or  predicted.  Davis 
and  McMorrow  (1976)11  and  Grantham  et  al.  (1983)12  provided  various  stochastic  models  for  deriving 
instantaneous  rainfall  from  clock-hourly  rainfall  for  selected  sites. 

From  available  observations,  one  may  express  rainfall  accumulation  in  terms  of  cumulative  percent 
frequency  per  rain  rate  for  any  unit  of  time.  Hydrologists  have  long  used  such  rainfall  mass  distributions 
for  estimating  maximum-depth  storms  or  flooding  possibilities  (World  Meteorological  Organization 
1973). 13 

Olascoago  (1950)14  and  Martin  (1964)15  found  that  when  daily  rainfall  per  year  is  converted 
to  cumulative  percent  frequency  per  cumulative  percent  amount,  despite  diverse  circumstances  the 
results  almost  invariably  tend  to  approach  a  particular  skew  distribution  (figure  1).  However,  Wexler 
(1967) 16  and  Cobb  (1968) 17  noted  that  this  average  or  nearly  universal  curve  does  not  describe  all 
rainfall,  particularly  extreme  climatic  conditions  or  rainfall  singularities.  Instead,  a  family  of  curves 
is  required  to  encompass  a  range  of  rainfall  accumulations,  whether  daily,  hourly,  or  instantaneous. 

The  main  purpose  of  this  paper  is  to  present  two  simple  models  for  rapid  estimates  of  short-term 
rainfall  for  any  unit  of  time.  These  are  a)  the  general  model  (cumulative  percent  frequency  per 
cumulative  percent  amount)  and  b)  the  explicit  model  (cumulative  percent  frequency  per  rain  rate). 
Each  model  consists  of  a  number  of  distributions,  which  depend  on  the  respective  average  rain  rate 
per  given  time. 

In  addition  to  a  viable  rainfall  mass  distribution  for  the  situation  at  hand,  the  essential  data  consist 
of  the  total  amount  of  precipitation  (P)  and  the  actual  duration  of  rain  in  days  (D),  hours  (H),  or 
minutes  (M)  for  estimating  daily,  hourly,  or  instantaneous  rainfall,  respectively.  A  climatic  index 
simulates  the  number  of  hours  of  rain,  if  necessary,  for  limited  regions.  The  only  instantaneous  rainfalls 
analyzed  are  for  single  storms  or  for  hourly  rain.  The  equations,  graphs,  and  computer  programs 
that  are  provided  facilitate  the  derivation  of  short-term  rainfall  for  a  wide  range  of  situations. 


'^Robert  O.  Berthel  and  Vernon  G.  Plank.  Time  Durations  of  Rain  Rates  Exceeding  Specified  Thresholds  (Hanscom  AFB.  Massachusetts: 
U.S.  Air  Force  Geophysics  laboratory,  1485,  4FGL-TR-0I22) 

*  *  Allen  R.  Davis  and  Daniel  J.  McMorrow.  Stochastic  Models  tor  Deriving  Instantaneous  Precipitation  Rate  Distributions.  Air  Weather 
Service  (Scott  AFB.  Illinois:  U.S.  Air  Force  Environmental  Technical  Applications  Center.  1976,  AWS  TR-76-263). 

*3D.D.  Grantham  et  al,"  Wafer  Vapor.  Precipitation.  Clouds  and  Fog:  Chapter  16.  1983  Revision,''  Handbook  of  Geophysics  and  Space 
Environments  (Hanscom  AFB.  Massachusetts  U.S  Air  Force  Geophysics  Laboratory.  July  1983.  AFGL-TR-83  0181). 

* 3 World  Meteorological  Organization.  "Manual  for  Estimation  of  Probable  Maximum  Precipitation.  "  Operational  Hydrology  ,  report 
no.  1  (Geneva:  1973.  WMO  no.  332). 

*4M.J.  Olascoago, ' 'Some  Asnects  of  A'gentma  Rainfall.  "Tellus  2  (November  1950)  3)2 

l5lA.  Martin.  "An  Investigation  of  the  Rainfall  Distribution  for  Stations  m  S'orth  and  Central  America.  "  Research  on  Tropical  Rainfall 
Patterns  and  Associated  Mesoscalc  Systems,  report  no  5  ((  ■liege  Station  Texas  A&M  University.  Dept  of  Meteorology,  contract 

DA-IAO-1 1001 -B-02 1-09.  May  1964) 

**Ruth  L.  Wexler,  "The  Relationship  fvrueen  Daily  and  Annual  Ramljil  in  Thailand.  Proceedings  of  the  Technical  Exchange  Conference. 
Air  Weather  Service  technical  report  196  (Washington.  DC  An  Wcathe:  Service.  C  S  Air  Force.  |467i  226  242. 

'7L.G  Cobb.  "The  Annual  and  Daily  Distribution  of  Rainfall  in  Southeast  Asia.  Resc.ir,  a  on  T  topical  Rainfall  Patterns  and  Associated 
Mesoscale  Systems,  report  no  4  (College  Station  Texas  A  AM  University.  Dept  .>!  Meteorology .  contract  DA  043-AMC  02313(Fl. 

December  1968)  53  75 


The  data  utilized  are  from  a  variety  of  sources,  including  National  Energy  Authority,  Thailand, 
1962-1965;  U.S.  Army  Engineer  Division  1968;  Sterntein  1962;  Environmental  Technical  Application 
Center  (ETAC)  1962;  U.S.  Weather  Bureau  1960;  Winner  1968;  and  Berthel  and  Plank  1985  (see 
references). 


THE  GENERAL  MODEL 


The  general  model  is  considered  first  because  of  its  overall  applicability.  The  example  of  the 
general  model  in  figure  1,  referred  to  as  the  ’’universal’’  or  default  rainfall  distribution,  is  used  to 
explain  the  methodology.  This  distribution  confirms  everyday  observations  that  rain  tends  to  fall 
unevenly,  with  heavy  rain  lasting  briefly  compared  to  gentler  rain.  More  specifically,  this  figure 
indicates  that  about  10  percent  of  the  rain  falls  in  about  50  percent  of  the  time,  whereas  nearly  50 
percent  of  the  rain  falls  in  just  about  10  percent  of  the  time. 

By  means  of  a  satisfactory  equation  for  the  plot  in  figure  1 ,  one  may  approximate  the  original 
short-term  rainfall  intensities  if  the  total  rainfall  and  the  duration  of  the  rain  are  known.  Unless  otherwise 
stated,  the  percent  frequency  refers  to  the  percent  of  the  actual  rain  days  (hours,  minutes).  One  possible 
equation  for  the  distribution  in  figure  1  is  of  the  form 

y  **  A  +  B(ln  x)  (1) 

where  y  =*  cumulative  percent  frequency 

x  =  cumulative  percent  amount  of  rain 
and  A  and  B  are  constants 

A  best-fit  equation  for  figure  1  is 

A  =  3.6 
B  =  21 


or  y  =  3.6  +  21(ln  x)  (la) 

dx/dy  — ^  =  x/B 

let  S  -  dx/dy  (lb) 

The  technique  is  expounded  in  terms  of  daily  rainfall,  but  the  method  is  applicable  to  any  short-term 
rainfall  (hourly  or  instantaneous)  as  shown  below.  Alternate  values  for  the  constants  A  and  B  are 
given  later  for  more  extreme  rainfall,  whether  very  light  or  very  heavy. 


3 


At  any  percent  frequency,  y,  the  corresponding  rain  rate,  r,  may  be  found  as  follows 


let 

then 

but 

therefore 

and 

where 


z  =  (y-A)/B 
x  =  exp(z) 

S  =  x/B 
r  =  S  (P/D) 
d  =  (100-y)/100 

P  =  total  rainfall  (mm  or  inches) 

D  =  number  of  days  with  measurable  rain 

r  =  rain  rate  mm  or  inches  per  day,  (compatible  with  P) 

d  =  the  number  of  days  with  rain  greater  than  r 


(lc) 

(ld) 


(For  hourly  rainfall,  substitute  H  and  h  for  D  and  d,  respectively.) 

Figure  2,  a  graphic  aid,  i.e.,  a  plot  of  S  versus  y,  provides  a  quick  way  to  estimate  short-term 
rainfall  with  essentially  little  computation.  For  instance,  for  any  y,  read  off  S,  then  multiply  by  (P/D) 
to  obtain  r.  As  shown  later,  equation  (la)  tends  to  underestimate  r  for  values  of  y  greater  than  95 
percent,  for  which  cases  extreme- value  statistics  are  indicated  (United  States  Weather  Bureau  1955; 
Hershfield  1972;  Fletcher  and  Sarlos  1951). 18  1920 


Inasmuch  as  P/D  (the  ratio  employed  in  most  precipitation  models)  provides  a  clue  to  the  nature 
of  short-term  rainfall,  a  question  raised  is,  How  frequently  is  this  average  rain  rate  (P/D)  equalled 
or  exceeded?  If  r  =  P/D  in  equation  (lc). 


Then  S  =  1 
x  =  B 

and  y  =  67.5% 

or  100-y  =  32.5% 


I&XJ.S.  Weather  Bureau.  "Run fill  Intensity- Duredon-Freqoency  Curves  for  Selected  Stetsons  in  the  United  Steles,  Aleske,  Hsweuen 
blende,  end  Puerto  Rico,"  technical  paper  no.  23  (1933). 

19D  M.  Henhfteld,  “ Esamedng  the  Extreme-  Value  1-Minute  Reinfell,  "  Journal  of  Applied  Meteorology  1 1  (1972):  936-940 
2®R.D.  Fletcher  and  D  Sarloa,  World  Record  Reinfell  end  en  Envelope  of  World  Record  Vetoes.  Air  Weather  Service  technical  report 
no.  103-SI  (Waihingtcn,  D  C  .  Air  Weather  Service,  I  S.  Air  Force,  1951). 
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If  P/H  were  substituted  for  P/D  in  the  above  equations,  the  answer  would  be  the  same  for  hourly 
rainfall.  Bussey  (1950)21  found  that  for  Washington,  D.C.,  the  average  hourly  rain  rate,  P/H,  is 
equalled  or  exceeded  about  35  percent  of  the  rain  hours. 

By  means  of  equations  (la)  to  (Id),  short-term  rainfall  is  estimated  for  various  situations  (table 
1).  Tables  la  and  lb  compare  observed  and  estimated  daily  rainfall  frequencies  per  year  and  month, 
respectively.  Table  lc  is  similar  for  hourly  frequencies  per  year.  The  results  appear  to  be  fairly 
reasonable,  considering  the  fact  that  all  the  above  estimates  are  made  by  means  of  the  distribution 
in  figure  1 .  Further  implications  of  the  results  are  discussed  in  more  detail  later. 

Equation  (la)  is  one  of  several  equations  that  describe  the  plot  in  figure  1 .  Another  is  the  exponential 

x  =  a  exp(by)  (2) 

where  a  and  b  are  constants 
a  «  1.03 
b  =  .0455 

x  and  y  are  the  same  as  above 

then  x  =  exp(.046y)  (2a) 

dx/dy  =  .046  x 

and  S  *  x/22  (2b) 

The  constant  a  is  rounded  off  to  a  value  of  1 .0  in  order  to  simplify  the  equation  (2a),  which,  therefore, 
approximates  but  does  not  exactly  duplicate  equation(la).  Appendix  A  contains  sample  problems 
that  are  solved  by  either  equation(la)  or  (2a). 

Despite  fairly  reasonable  estimates  of  daily  or  hourly  rainfall  for  the  temperate  zones  and  the 
tropics  by  means  of  the  universal  curve  (figure  1),  this  distribution  has  its  limitations,  as  mentioned 
above,  particularly  with  respect  to  relatively  light  rain,  characteristic  of  the  very  high  latitudes,  or 
very  intense  rain,  as  at  certain  desert  stations,  or  in  very  brief  storms.  Figure  3,  therefore,  consists 
of  a  family  of  curves  (in  the  format  of  figure  1)  encompassing  various  rainfall  occurrences,  as  extreme 
climates  or  intense  storms.  These  curves  represent  daily  and  hourly  as  well  as  instantaneous  rainfall 
frequencies  based  on  more  than  100  cases.  Each  curve  is  a  smoothed  average  of  several  cases  except 
for  the  dotted  curves,  which  are  only  for  a  single  case  or  two.  In  a  given  situation,  a  particular 
distribution  might  approach,  but  not  actually  duplicate,  any  of  the  above  curves.  More  than  half  the 
cases  analyzed  consist  of  original  observations  of  daily  or  hourly  rainfall  per  year  or  per  month. 


“'Hn.vaid  F.  Bussey.  "M/crowive  Attenuation  Statistics  Estimated  from  Rainfall  and  Water  Vanor  Statistics. "  Proceedings  of  the  Institute 
of  Radio  Engineers  38  (New  York:  1950):  781-785. 
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The  rest  are  of  tabulated  class  frequencies  or  previously  plotted  data  of  daily  or  hourly  rainfall  per 
year  or  instantaneous  rainfall  per  single  storm.  Curve  No.  1 ,  the  leftmost  curve  in  figure  3,  represents 
hourly  rainfall  for  a  single  year  at  Pensacola,  Florida.  The  majority  of  cases  fall  within  the  envelope 
formed  by  curves  No.2  and  No.4,  with  curve  No.3  (bold  line)  being  essentially  the  same  as  in  figure 
1.  The  more  extreme  the  rainfall,  the  greater  the  departure  from  the  latter  distribution.  Table  B1 
(appendix  B)  contains  constants  for  equation  (1)  for  the  different  curves  in  figure  3.  If  rainfall 
accumulations  are  identified  by  their  respective  average  rates,  as  P/D,  P/H,  or  P/M,  some  tentative 
distributions  based  on  the  numbered  curves  in  figure  3,  are  as  follows: 


Curve 

Period 

2-4 

year 

year 

5 

year 

year 

1  hour 

6 

1  hour 

Average  Rain  Rate 
P/D:  8.2  to  21.5  mm/day 
P/H:  1.2  to  4.0  mm/hour 
P/D  >  28.0  mm/day 
P/H  >  25.0  mm/day 
P/M  >  50.8  mm/hour  (0.9  mm/minute) 
P/M  >  127  mm/hour  (2.1  mm/minute) 


Curves  No. 5  and  No. 6  both,  then,  represent  (approximately)  the  distributions  of  selected  maximum- 
depth  one-hour  storms,  No.5  for  50  mm  (2  inches)  and  No. 6  for  127mm  (5  inches).  For  the  first 
storm,  30  percent  of  the  rain  occurs  in  10  percent  of  the  time  (6  minutes),  whereas  for  the  second 
storm,  22  percent  of  the  rain  occurs  in  10  percent  of  the  time.  If  these  6-minute  intervals  were 
converted  to  hourly  rates,  the  results  would  be  152  mm  and  279  mm,  respectively.  The  latter  rate 
is  about  three-fourths  of  world-record  rainfall  for  1  hour  (381  mm),  as  found  by  Fletcher  (1951). 22 
For  values  of  y  equal  to  or  greater  than  70  percent  or  for  distributions  with  curve  numbers  greater 
than  5,  an  alternate  equation  is 


y  =  a  xb  (3) 

Appendix  B  contains  the  respective  constants  for  equation  (3)  for  the  above  two  maximum- 
depth  one-hour  storms. 


22r  d.  Fletcher  and  D.  Sarlos,  World  Record  Rainfall  and  an  Envelop  of  World  Record  Values  Air  Weather  Service  technical  report 
no.  105-8 1  (Waahingtnn.  D.C.,  Air  Weather  Service,  U.S.  Air  Force,  1951). 
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THE  EXPLICIT  MODEL 


The  explicit  model  consists  of  a  series  of  rainfall  distributions  expressed  in  terms  of  cumulative 
percent  frequency  per  rain  rate,  as  in  figures  4  and  S.  Unless  otherwise  indicated,  each  distribution 
is  the  average  of  3  or  more  years  of  annual  rainfall  data  for  a  given  station.  In  each  figure,  the 
selected  stations  represent  a  range  of  daily  or  hourly  rainfall  regimes.  Although  the  selection  was 
made  at  first  on  the  basis  of  a  wide  spectrum  of  total  rainfall  per  year,  the  resultant  plots  are  not 
necessarily  ordored  by  rainfall  amount. 

Figure  4  compares  daily  (calendar-day)  rainfall  distributions  among  stations  in  Thailand,  with 
total  rain  from  about  1000  mm  to  over  3000  mm  per  year.  An  interesting  feature  of  this  figure 
is  the  increasing  value  of  P/D  from  left  to  right  regardless  of  the  individual  values  of  P  or  D. 
For  P/D  of  8  mm/day,  less  than  10  percent  of  the  rain  days  are  more  than  23  mm,  whereas  for 
P/D  of  31.8  mm/day,  nearly  30  percent  of  the  rain  days  are  more  than  23  mm.  Table  2  lists  the 
original  stations  and  the  corresponding  annual  rainfall  represented  by  these  plots. 

Figure  3,  for  hourly  (clock-hour)  rainfall,  shows  dramatically  the  variation  of  rainfall 
distributions  from  the  Arctic  to  the  Tropics.  As  with  daily  rainfall,  the  average  rain  rate  in  figure 
3,  i.e.  P/H,  increases  from  left  to  right  irrespective  of  P  or  H.  In  the  high  latitudes  (P/H  =  ,6mm/hr), 
just  a  small  percentage  of  the  rain  is  more  than  2.3  mm/hr  (.  1  in/hr),  whereas  in  the  low  latitudes 
(P/H  *  4.0  mm/hr  to  nearly  6  mm/hr),  40  percent  of  the  rain  is  more  than  2.3  mm/hr.  At  Roi 
Et,  Thailand,  and  Veunesai,  Cambodia  (not  shown  in  figure  3),  both  with  P/H  =  5.0  mm/hr, 
about  30  percent  of  the  rain  hours  are  more  than  2.3  mm,  and  about  10  percent  of  the  hours  are 
more  than  12.7  mm.  Table  3  lists  the  stations  and  corresponding  annual  hourly  rainfall  data  for 
die  plots  in  figure  4. 

Table  4  lists  various  percentiles  of  instantaneous  rainfall  per  hour  for  essentially  continuous 
rain  for  average  hourly  rates  of  from  0.42  mm/min  to  2.12  mm/ min  (1  inch  to  5  inches  of  rain). 
These  data  are  from  records  of  excessive  rain  for  periods  of  3  minutes  to  60  minutes  (Panama 
Canal  Companyl957-1971).23 

One  equation  for  the  distributions  in  figures  4  and  5  or  table  4  is  similar  to  equation  ( 1 )  except 
that  R  replaces  x,  that  is. 


y  -  A  +  B(ln  R)  (4) 

where  R  *  rain  rate 

Inasmuch  as  each  of  the  constants  A  and  B  appears  to  be  a  function  of  P/D,  P/H,  or  P/M, 
intermediary  values  of  either  constant  may  be  interpolated  accordingly  (appendixes  B  and  C).  Tables 
2  to  4  contain  values  of  the  constants,  A  and  B,  for  equation  (4)  for  the  distributions  indicated. 


23Panama  Canal  Company,  Climatological  Data,  Canal  Zone  and  Panama,  annual  reports  1957-1971  (Balboa  Heights,  Canal  Zone). 
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As  in  the  case  of  the  general  model,  the  cumulative  percent  frequency,  y,  of  the  average  rain 
rate,  P/D  (P/H),  which  identifies  the  distribution,  may  be  readily  computed  by  means  of  the  explicit 
model,  this  time  by  substituting  P/D  (P/H)  for  R  in  equation  (4).  For  example,  for  P/D  =  21.5 
mm/day,  A  =  -27,  B  =  30.0,  and  y  =  65%  (see  figure  4). 

A  single  distribution  may  not  be  as  reliable  for  the  explicit  model  as  a  series  of  distributions 
for  a  range  of  P/D  or  P/H  (or  P/M )  from  which  the  appropriate  constants  A  and  B  may  be 
interpolated  if  necessary.  The  term”explicit”  implies  explicit  per  given  P/D  (P/H  or  P/M)  in  contrast 
to  the  general  model,  where  a  single  distribution  may  serve  a  considerable  range  of  daily,  hourly, 
or  instantaneous  rainfall. 

Since  the  nature  of  the  distribution  of  short-term  rainfall  per  given  situation  appears  to  depend 
on  the  respective  average  rain  rate  (P/D,  P/H,  or  P/M)  for  the  interval  in  question,  some  previous 
observations  are  now  reexamined  (tables  5  to  7).  For  instance,  table  5  (including  some  of  the  annual 
daily  rainfall  data  of  table  1)  lists  43  Thai  stations,  with  P/D  (per  year)  varying  from  8.9  to  24.5 
mm/day,  respectively,  from  10-year-average  data.  The  original  observations  (from  Sternstein 
1962)24  have  been  converted  for  this  purpose  to  cumulative  percent  frequency  per  rain  rate.  Table 
6  lists  daily  and  hourly  frequencies  per  season  and  year:  6a  is  in  conventional  format;  6b  and 
6c  are  in  terms  of  percent  frequencies  of  the  actual  rain  duration,  by  order  of  P/D  and  P/H, 
respectively.  The  latter  tables  show  that  the  cumulative  percent  frequency  per  rain  rate  changes 
very  slowly  with  increasing  P/D(P/H).  This  fact  no  doubt  accounts  for  the  apparent  utility  of  the 
universal  (default)  model.  Each  of  the  above  tables  serves  as  an  explicit  model.  Table  6c  indicates 
that  seasonal  rainfall  frequencies  follow  the  annual  trends,  but  have  greater  disparity  in  their 
distributions.  Moreover,  this  table  shows  that  two  cases  of  seasonal  or  annual  rainfall  with  the 
same  P/H  do  not  necessarily  have  the  same  P/D  (6b),  or  vice  versa.  For  instance,  San  Juan-Summer 
and  Hilo-Winter  both  have  a  P/H  of  2.26  mm/hr,  whereas  the  respective  P/D’s  are  8.31  and  15.42 
mm/day.  The  relationship  between  H  and  D(or  H  and  P)  apparently  has  to  be  established  for  the 
situation  at  hand. 

For  each  unit  of  time,  as  a  day  or  an  hour,  there  appears  to  be  a  limiting  rain  rate  beyond 
which  the  cumulative  percent  frequencies  scarcely  vary  with  increasing  P/D  or  P/H.  For  example, 
in  table  5,  from  P/D  of  11.4  mm/day  to  14.9  mm/day,  the  frequency  for  rain  of  35  mm  a  day 
is  about  90  percent.  For  P/D  of  8.9  mm/day  to  14.9  mm/day,  for  a  daily  rain  of  90  mm,  the 
frequencies  are  all  more  than  99  percent.  In  table  7,  the  frequencies  associated  with  hourly  rain 
greater  than  12.7  mm  (.50  inch)  are  similar  for  values  of  P/H  from  2.8  mm/hr  to  5.0  mm/hr. 
As  a  consequence,  the  following  formula  is  utilized  for  given  limits  of  y  or  R,  namely, 

y  =  a  Rb  (5) 

See  appendix  B  for  constants  for  equation  (5). 


Stenutein.TTie  Rtinfal]  of  Thailand  (Bloomington:  Indiana  University  Foundation  Research  Division,  contract  DA-I9-129-QM-1582, 

1962). 
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UTILITY  OF  THE  DIFFERENT  MODELS 


In  addition  to  tables  la  to  lc  and  appendix  A,  which  give  examples  of  the  use  of  the  universal 
distribution  in  figure  1,  appendix  C  lists  a  number  of  computer  programs  that  show  a  variety  of 
ways  for  which  both  models  may  be  exploited.  Unless  otherwise  stated,  only  equations  (la)  or 
(2a)  are  used  to  demonstrate  the  general  model  (the  exception  being  instantaneous  rainfall,  or  the 
last  2  programs). 

The  computer  programs  listed  below  determine  the  following,  per  P  and  D,  H,  or  M,  for 
the  case  at  hand,  namely, 


1)  DAYRATE:  Estimates  die  number  of  days  per  year  (month,  season)  that  a  given  rain  rate  is  equalled 
or  exceeded. 

2)  HOURRATE:  Estimates  the  number  of  hours  per  year  (month,  season)  that  selected  rain  rates 
are  equalled  or  exceeded. 

3)  FREQHOUR:  Determines  the  rain  rate  that  is  equalled  or  exceeded  0.1,  0.2,  0.S,  1.0,  and  2.0 
percent  hours  per  year. 

4)  TESTDAY :  Compares  the  explicit  and  the  general  models  for  daily  rain. 

3)  TESTHOUR:  Compares  the  explicit  and  the  general  models  for  hourly  rain. 

6)  INSTANT:  Determines  the  deciles  of  instantaneous  rainfall  for  selected  maximum-depth  one-hour 
storms  in  terms  of  mm/min. 

7)  ONE-HOUR:  Computes  the  distributions  of  instantaneous  rainfall  for  several  one-hour  storms 
in  terms  of  mm/hr. 


AN  HOURLY  RAINFALL  CLIMATIC  INDEX 


Inasmuch  as  the  number  of  hours  of  rain  (H)  is  not  usually  available,  it  must  often  be  obtained 
by  indirect  means.  Table  8  lists  P,  D,  and  H  for  selected  stations.  On  a  worldwide  basis,  H/D  * 
4.0,  approximately.  However, this  ratio,  somewhat  erratic,  tends  to  be  higher  in  the  colder  regions 
than  in  the  warmer.  From  the  high  to  the  low  latitudes,  P/H  varies  from  0.6  mm/hr(.02  inches/hr) 
to  about  10  times  that  rate,  respectively.  Let  J  =  H/P:  then  for  any  annual  rainfall,  J  may  be  roughly 
estimated  (t2=.16)  from  mean  annual  temperature  (°C)  by 
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J  =  .0447(C)  +  1-43 


(6) 


where  J  is  in  hours/mm 

Equation  (6)  does  not  apply  to  stations  with  P/H  greater  than  4.0  mm/hr.  Once  J  is  established  for 
a  particular  climatic  zone,  then  for  any  station  in  the  respective  area,  H  =  P(J). 

In  any  given  region,  the  value  of  J  may  be  locally  influenced  by  exposure  to  the  moist  air  flow 
or  by  terrain.  For  rainfall  in  Indonesia  (Mohr  and  VanBaren  1954)25,  for  instance,  J  tends  to  increase 
somewhat  with  elevation. 


DISCUSSION 


The  greatest  advantage  of  the  general  model  is  its  wide  applicability,  that  is,  the  use  of  essentially 
a  single  distribution  for  a  considerable  variety  of  circumstances.  Moreover,  for  a  given  situation, 
the  constants  for  the  general  model  remain  the  same  no  matter  what  the  data  units  (inches/day  or 
mm/day).  For  the  explicit  model,  the  constant  A  has  to  be  converted  with  change  of  unit,  whereas 
B  remains  the  same.  For  instance,  A(inches)=A(mm)  +B(ln  25.4).  For  either  model,  the  relative 
stability  of  rainfall  distributions  in  terms  of  percentiles  provides  a  means  of  checking  data  errors. 
The  main  disadvantage  of  the  general  model  is  that  unless  the  data  represent  a  long-term  average, 
a  single  extreme  rain  rate  on  occasion  may  considerably  bias  the  resultant  distribution. 

The  advantage  of  the  explicit  model  is  that  one  may  sometimes  estimate  rain  rate  frequencies 
merely  by  inspection  (or  interpolation)  of  the  appropriate  graph  or  table.  This  model  provides  a  guide 
to  the  limiting  values  of  the  short-term  rainfall  intensities  (per  rainy  period)  in  the  natural  environment 
(figures  4  and  5).  The  drier  the  climate  (or  season),  the  greater  the  dearth  of  low  intensity  daily 
(or  hourly)  rainfall.  For  instance,  for  3  years,  1962-1964,  of  rainfall  at  Sisaket,  Thailand,  with  P/D 
values  of  38.4,  27.5,  and  28.2  mm/day  and  D  of  50,  54,  and  30,  respectively,  there  was  not  a  single 
rain  day  with  less  than  4  millimeters  of  rain.  The  colder  the  climate  (or  season),  the  greater  the 
percentage  of  very  light  rainfall  per  day  or  hour.  For  annual  daily  rainfall,  a  limiting  distribution 
(see  figure  4)  is  reached  at  about  30  mm/day  (calendar-day).  Any  higher  average  daily  rain  rate  often 
(not  always)  results  from  a  single  day  of  highly  unusual  rainfall.  For  annual  clock-hourly  rainfall, 
the  limit  is  about  5  mm  per  hour  (see  table  7).  For  continuous  (instantaneous)  rainfall,  some  measure 
of  the  total  range  of  the  instantaneous  intensities  seems  to  be  required  (Berthel  and  Plank  1985)26, 
a  factor  which  might  also  improve  extreme  daily  or  hourly  rainfall.  For  one-hour  storms  of  50.8 
mm  and  127  mm,  the  range  of  instantaneous  rates  appears  to  be  greater  in  the  U.S.  Weather  Bureau 
examples  than  at  the  Canal  Zone  and  Panama  (compare  the  results  of  the  programs,  INSTANT  and 
ONE-HOUR,  appendix  C). 


t  Mohr  and  F  A.  Van  Barer,  Tropical  Soils  (London:  Interscience  Publishers,  Lid.,  1954) 

"Robert  O.  Berthel  and  Vernon  G.  Plank,  Time  Durations  of  Rain  Rates  Exceeding  Sepcified  Thresholds  (Han  scorn  AFB,  Massachusetts: 
U  S.  Air  Force  Geophysics  Laboratory,  1985,  AFGL-TR-0122) 
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SUMMARY  AND  CONCLUSIONS 


1 .  Despite  the  vagaries  of  rainfall  and  the  diversity  of  circumstances,  for  any  one  year  or  for  a  single 
storm  die  rainfall  distributions  determined  so  far  are  surprisingly  reproducible  in  terms  of  cumulative 
percent  frequency  per  cumulative  percent  amount  or  else  cumulative  percent  frequency  per  rain  rate, 
for  a  given  P/D  (P/H  or  P/M),  regardless  of  the  individual  values  of  P  and  D  (H  or  M). 

2.  As  a  consequence  of  die  above  observations,  equations,  graphs,  and  computer  programs  are  presented 
to  facilitate  the  recovery  of  a  variety  of  information  concerning  short-term  rainfall,  particularly  its 
frequency. 

3.  Although  the  ratio  P/D  (P/H  or  P/M)  is  the  most  significant  factor  in  the  determination  of  short¬ 
term  rainfall  frequencies,  the  actual  duration  of  rain,  that  is  D  (H  or  M),  is  essential  for  converting 
the  percent  frequencies  to  durations  in  real  time. 

4.  The  hourly  rainfall  climatic  index,  H/P  or  J,  seems  to  be  a  useful  tool  for  deriving  H  when  only 
P  is  known. 

3.  For  annual  rainfall,  the  average  rain  rate,  P/D  (or  P/H),  is  exceeded  about  32  percent  to  35  percent 
of  the  rain  period. 

6.  Short-term  rainfall  distributions  tend  to  become  less  skewed  as  the  average  rain  rate  per  given 
interval  of  time  approaches  extreme  values. 

7.  In  general,  hourly  rainfall  intensities  increase  with  temperature  and  decrease  with  elevation. 

8.  Since  short-term  rainfall  has  implication  for  soil  moisture,  field  operations,  electro-optical  sensors, 
and  equipment  design  or  malfunction,  die  above  results  should  provide  invaluable  engineering  guidance. 
Military  planners  should  benefit  by  improved  understanding  of  the  likely  impact  of  rainfall  over  short 
periods  of  time  with  respect  to  mobility  or  trafficability. 


Cumulative 


Note:  For  the  percentile  of  any  daily  rain,  read  off 
S  at  selected  y,  then  multiply  S  by  P/D  to  obtain  the 
required  rate,  r.  Similarly,  for  hourly  rain,  multiply 
S  by  P/H. 
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Figure  3.  Set  of  rainfall  distributions. 
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5.  Cumulative  percent  frequency  of  hourly  rain 


DAILY  RAIN  FREQUENCIES  PER  YEAR:  THAILAND 
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short-term  rainfall 


DAILY  RAINFALL  FREQUENCIES  PER  MONTH 
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Table  1.  Observed  and  estimated  short-term  rainfall,  (cont) 
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Table  1.  Observed  and  estimated  sbort-tenn  rainfall,  (con) 
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Ikble  2.  Dally  rainfall  frequency  observations. 
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*Constants  for  equation  (4) 
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OBSERVATIONS  IN  ORIGINAL  FORMAT 
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Tkble  6.  Seasonal  daily  and  hourly  rainfall:  5  stations. 
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Table  6.  Seasonal  daily  and  hourly  rainfall:  5  stations,  (cont) 


CUMULATIVE  PERCENT  FREQUENCY  PER  HOURLY  RAIN  PER  P/H 
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Thble  6.  Seasonal  daily  and  hourly  rainfall:  5  stations,  (cont) 
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Table  8.  Daily  and  hourly  annual  precipitation  data  per  latitude. 
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APPENDIX  A 


Sample  Problems 


This  section  provides  sample  problems  and  their  solutions  by  means  of  equations  (la)  to  ( lc) 
and/or  figures  1  and  2. 

1)  Objective:  to  find  the  frequency  in  days  (d)  per  year  with  rain'rate,  r, 
>35  mm/day. 

Station  =  Chanthaburi,  Thailand 


P  =  3305  mm  (annual  rainfall) 
D  =  171  days  of  rain  per  year 
r  *  35 

S  =  35/(P/D)  -  1.81 
x  =  21.0(S)  *  38.0 
y  -  3  6  +  21(ln  38)  *  80.0 
d  =  D(100-y>/100  =  34 


Alternate  method  by  figure  2 

For  S  -  1.81,  read  (100-y)  =  20 
d  *  171(20/100) 


34  days  Ans. 
(obscrvea  =  29  days) 


or  d  =  34 


2)  Objective:  to  find  the  percent  hours  per  year  that  a  rain  rate  of  0.1  inch  per  hour  is  equalled 
or  exceeded. 

Station  =  Ceylon,  India 

P  =  90  inches  (annual  rainfall) 

H  =  607  hours  per  year  of  rain 
F  =  8760  total  hours  in  a  year 
r  =  0.1  inch/hr 
S  =  r/(P/H)  =  .67 
x  =  21(S)  =  14.1 
y  =  3.6  +  21(ln  x)  =  59 
h  =  H(100-y)/100  =  249  hours 

f  =  100(h/F)  =  2.8  2.8%  A  ns. 

(observed  =  2.7%) 

3)  Objective:  to  find  a  rain  rate,  r,  that  is  equalled  or  exceeded  2%  hours  per  year. 

Station  =  Washington,  Dt  C 

P  =  41.9  inches  (annual  rainfall) 

H  =  635  hours  of  rain  per  year 
h  =  175.2  =  2%  hours  per  year 
y'=  100(h/H)  =  27.6%  of  the  rain  hours 
y  =  100-y '  =  72.4 
x  =  26  (read  from  figure  1) 

S  =  x/21  =  1.24 

r  =  S(P/H)  =  .084  inch/hr  084  inch/hr  Ans. 

(observed  =  .076  inches/hr) 


4)  Objective:  to  find  the  percent  frequency  (f)  of  days  (d)  per  year  with  a  rain  rate  between  35 
and  10  mm/day. 

Station  =  Mae  Hong  Son,  Thailand 
P  =  1255  mm  (annual  rainfall) 

D  =  112  days  per  year 
r'  =  35 
r  =  10 

S'  =  r^P/D)  =  3.12 

y'  =  90  (read  y  at  S  =  3.12  on  figure  2) 

S  *  ^(P/D)  =  .89 

y  =  64.5  (read  y  at  S  =  .89  on  figure  2) 
d  =  (y'-yXD/100)  «  28.6  days 
f  *  100(28.6/365)  =  7.82%  days/year 

7.82%  Ans. 
(observed  =  8.5%) 


5)  Objective:  To  estimate  maximum  daily  rainfall 
Station  -  Sakon  Nakhon,  Thailand  (Sept  1964) 

P  =  262.7  mm  (monthly  rainfall) 

D  =  20  days  per  month 
f  =  100/D  =  5%  (1  day  of  rain) 
y  =  95 

S  =  3.7  (at  y  =  95  on  figure  2) 
r  =  S(P/D)  =  48.6 

48.6  mm  Ans. 
(observed  =  47.9  mm) 
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APPENDIX  B 


Additional/AItemate  Formulas 


1)  Explanation 

Since  neither  model  is  a  perfect  fit  and  since  rainfall  never  reproduces  itself  exactly,  the  equations 
given  at  best  provide  a  guide  to  common  trends  noted  in  long-term  or  even  single  sets  of  observations 
of  daily,  hourly  or  instantaneous  rainfall,  respectively.  Usually,  no  single  equation  is  applicable 
to  the  entire  distribution,  but  for  the  most  part  either  the  general  or  explicit  model  will  recover 
rain  rates  up  to  about  the  95th  frequency  percentile. 

To  the  extent  possible,  the  equations  employed  have  a  coefficient  of  determination  of  r2  >  .90 
(r2  =  1.0  represents  a  perfect  fit,  r2  =  0  is  no  fit).  All  equations  are  subject  to  revision  as  data 
improve. 

Discrepancies  in  the  constants  come  about  if,  for  any  given  formula  and  data  set,  a  difference 
obtains  in  the  frequency  range  (y)  selected. The  observations  themselves,  as  figures  4  and  5  or 
tables  4,  5,  6b,  6c,  or  7b  provide  the  most  reliable  information  concerning  the  frequency  of  short¬ 
term  rainfall.  The  formulas  are  sensitive  to  the  size  of  the  class  frequencies  utilized. 

2)  Equations  for  x  vs  y 

where  y  =  cumulative  percent  frequency 
x  =  cumulative  percent  amount 

and  A,B,a,b  are  constants 

a)  Equation  (1):  y  =  A  +  B(ln  x) 

Table  B1  lists  constants  A  and  B  for  the  curves  in  figure  3  of  the  general  model. 

Of  the  distributions  in  figure  3,  curve  No.  3  (same  as  figure  1)  is  referred  to  as  the  universal 
or  default  distribution. 

For  instantaneous  rainfall,  from  a  set  of  depth-duration  plots  for  one-hour  storms  of  50  mm 
to  127mm  total  amounts  (U.S.  Weather  Bureau),  equations  of  the  above  format  have  been 
determined,  namely, 


50  mm: 

y  =  -35.1  +  28.9(ln  x) 

(Bl) 

127  mm: 

y  =  -69.7  +  35.8(ln  x) 

(B2) 
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Once  a  set  of  equations,  as  (Bl)  and  (B2)  has  been  obtained  for  two  or  more  cases  of  specified 
hourly  rainfall,  then  the  constants,  A  or  B,  for  any  intermediate  rainfall  may  be  readily  derived 
as  functions  of  P/H  or  P/M.  For  P/M,  from  0.83  mm/ min  to  2.12  mm/ min  (corresponding  to 
30  mm/hr  and  127  mm/hr,  respectively),  A  and  B  may  be  found  as  follows 


A  -  -12  -27.2(P/M)  (B3) 

B  *  24.3  +  5.5(P/M)  (B4) 

P/M  *  mm/min 

b)  For  equation  (2):  x  =  a  exp  (by) 

50  mm:  x  *  3.42  exp(.03  y)  (B5) 

127  nun:  x  =  7.14  exp(.03  y)  (B6) 

c)  For  equation  (3):  y  ■  a  xb 

50  mm:  y  »  7.51  x  “  (B7) 

127  mm:  y  -  3.57  x  *  (B8) 

3)  Equations  for  y  vs  R 

where  y  =  the  same  as  above 


R  »  rain  rale  (in  specified  units) 
a)  For  equation  (4):  y  =  A  +  B(ln  R) 

For  daily  rain,  8.2  mm/day  to  17.0  mm /day,  the  derivation  formulas  for  die  constants  A  and  B  are 

A  =  38.8-1. 8(P/D)  (B9) 

B  =  21  (BIO) 

For  hourly  rain,  2.8  mra/hr  to  4.5  nam/hr  (data  from  various  sources),  the  derivation  formulas 
for  the  constants  A  and  B  ate 


A  =  64.7  -6.6(P/H)  (r  =  .8) 
B  -  21.7 

b)  For  equation  (5):  y  =  a  Rb 

a  b 

For  y  >70  56.4  0.20 

For  R  values  (mm/hr) 

12.7  -  25.4  74.4  0.09 


(BID 

(B12) 
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■«t  *o  t- 


CURVE#  A  B 


30.7 

16.5 

16.1 

19.0 

3.6 

21.0 

-15.8 

25.0 

-37.2 

29.7 

-69.0 

36.2 

-107.0 

44.0 

Table  Bl.  Constants  for  equation  (1)  for  curves  in  figure  3. 
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APPENDIX  C 


Cwnpoter  Programs 


The  distribution  in  figure  1  (general  format)  is  utilized  as  the  default  or  universal  model  for  the 
first  five  programs.  However,  programs  No.4  and  No.S  also  include  the  explicit  model,  which  employs 
equations  (4)  and  (S).  The  required  constants,  A  and  B,  may  be  assigned  if  they  are  available  for 
the  rainfall  regime  in  question,  that  is,  the  particular  P/D,  or  P/H.  Otherwise,  they  may  be  derived, 
as  explained  previously. 

If  any  rainfall  regime  being  investigated  has  a  value  of  P/D  or  P/H  within  the  limits  of  those 
given  in  figures  4  or  5  or  tables  4  to  7,  very  often  reasonable  estimates  in  terms  of  percentiles  may 
be  made  by  inspection  for  comparison  with  computer  results. 
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1)  DAYRATE 


1  REM  "DAYRATE"  RLW  86 

2  LPRINT 

3  PRINT’  PROGRAM  TO  EST  DAILY  RAIN  RATE  PROBABILITIES" 

4  REM  EQUATION:  Y  -  3.6  ♦  21  LOG  (X)  NOTE  LOG  -  NAT  LOG 

6  REM  DX/DY  -  X/21  LET  S  -  DX/DY 

7  REM  THE  RAIN  RATE,  R,  AT  ANY  Y  IS  R-S*(P/D) 

8  LPRINT’STATION" ."INTERVAL" P(MM)  D  ";"  RATE(R)";"  DAYS>R 

9  LPRINT 

10  PRINT’ENTER  STATION" 

12  INPUT  S$ 

20  PRINT  "YEAR,  SEASON,  OR  MONTH?  IF  SEASON,  WHICH?  IF  MONTH,  WHICH?” 
22  INPUT  1$ 

30  PRINT  "ENTER  TOTAL  RAIN  IN  MILLIMETERS" 

35  INPUT  P 

70  PRINT  "  NUMBER  OF  RAIN  DAYS" 

80  INPUT  D 
82  A  -  3.6 

84  B  -  21 

85  I  -  P/D 

90  PRINT  "ENTER  RAIN  RATE  IN  MM" 

100  INPUT  R 

130  S  -  R/I 

140  X  -  S  *  B 

150  Y  -  A  +  ( B*LOG (X ) ) 

152  IF  Y <0  THEN  Y  -  OsGOTO  160 
154  IF  Y>1 00  THEN  Y  -  100 
160  F  ~(100-Y)*D/100 
164  IF  Y>1 00  THEN  F  -  0 

170  LPRINT  S$, 

171  LPRINT  1$, 

172  LPRINT  USING"####.#  "jP.D.R.F 

180  PRINT" ANOTHER  CASE?  Y?  N?" 

182  INPUT  A$ 

184  IF  A$»"Y"  THEN  GOTO  10 
190  END 


SAMPLE  OUTPUT:  DAYRATE 


STATION 

INTERVAL 

P  (  MM) 

D 

RATE(R) 

DAYS>R 

WASH . D . C . 

AUTUMN 

203.2 

24.0 

6.4 

9.2 

HILO 

SPRING 

617.2 

70.8 

12.7 

17.4 

BERLIN 

YEAR 

584.2 

171.0 

1.0 

99.6 

CHANTHABURI 

YEAR 

3304.5 

170.7 

10.0 

79.1 
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2)  BOOBEATB 

2  KKM  "BOOBBATS*  IB/BB  BLV  86 

4  IBM  TO  DBT  BO.  BOOBS  HXTB  BAXN>-.05,  .10,  .25  INCHES 
6  OPTION  BASS  X 

5  DIM  1(3) ,  B(3) 

10  LPBXMT 

11  LPBXMT 

12  LPBIMT" ”TAB(30)*M0 .  BOOBS  >  XN/BB  GIVEN" 

14  LPBIMT 

16  B(l)  -.05 

18  B(2)  -9.999999B-02 

20  B(3)  -.25 

22  LPBIMT* STATION1’,*  P(XM)  *;  *  B  *;*  .05 

23  LPBIMT*  .10  *;*  .25  * 

24  LPBIMT 

30  PBIMT*BMTBB  STATION* 

32  IMPOT  S$ 

34  PBIMT  "BMTBB  TOTAL  BAXM  XM  INCHES" 

36  IMPOT  P 

38  PBIMT  "BMTBB  TOTAL  BOOBS  OP  BAXM* 

40  IMPOT  M 
42  X  -  P/B 
44  FOB  L  -  1  TO  3 
50  S  -  B(L)/1 
52  X  -  S*21 

54  I  -  3.6+21*<L0C(X)) 

60  IP  (Y<0)  TMBM  T-0  t  GOTO  10 
62  IP(T>100)  TBBM  T-100 
70  P(L)  -  <100-T)«B/100 
72  MBIT  L 
80  LPBIMT  SS, 

90  LPBIMT  OSING*####.#  "jP,  B,  P(1 ) ,P(2 ) ,P(3 ) 

100  PBIMT" AMOTBBB  CASBf * 

110  IMPOT  AS 

120  IP  A|-*T*  TMBM  GOTO  30 
130  BMD 


SAMPLB  OOTPOTl  BOOBEATB 


MO.  BOOBS  >  XM/HB  GIVEN 


STATION 

KIM) 

8 

.05 

.10 

.25 

T0NMSV1LLB 

45.7 

364.0 

188.5 

135.6 

65.5 

PBSBAWAB 

14.7 

109.0 

58.1 

42.2 

21.3 

WASHINGTON 

41.9 

635.0 

243.1 

150.7 

28.5 

EBAKOE 

50.9 

388.0 

204.6 

148.1 

73.4 
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3)  FREQBOUB 


2  REM  "FREQBOUB  "  RLW  86 

3  REN  THIS  IS  AN  INTERACTIVE  PROGRAM 

4  REM  TO  FIND  RAIN  RATE  EQUALLED  OR  EXCEEDED 
6  REM  THE  Z  HOURS/TEAR  INDICATED 

8  REM  TOTAL  HOURS  PER  TEAR  -  8760;  1Z  OF  TIME  -  87.6  HOURS 
10  REM  PARAMETERS  REQUIRED:  P  -  TOTAL  ANNUAL  RAIN,  H  -  HOURS  OF  RAIN 
12  DATA  0.1,  0.2,  0.5,  1.0,  2.0 
14  DIM  R<5),  F(5 ) ,  C(5),  Y(5) 

16  FOR  L  -  1  TO  5 
18  READ  C(L) 

20  F(L )  -  C(L )  *  87.6 
22  NEXT  L 
24  LPRINT 

26  LPRINT"  "TAB( 38 ) "RAIN  RATE (MM/ HR)  EQUALLED  OR  EXCEEDED" 

27  LPRINT"  "TAB( 38 ) "  Z  HOURS  INDICATED" 

28  LPRINT 

30  LPRINT" STATION" ,"  P  H  "; 

32  LPRINT"  .1Z  ";"  .2Z  ";"  .5Z  1Z  2Z  " 

34  LPRINT 
40  A  -  3.6 
42  B  -  21 

44  PRINT" ENTER  STATION" 

46  INPUT  S$ 

SO  PRINT" ENTER  TOTAL  ANNUAL  RAIN  IN  MILLIMETERS" 

52  INPUT  P 

54  PRINT" ENTER  NUMBER  OF  HOURS  OF  RAIN  PER  YEAR" 

56  INPUT  H 
60  FOR  L  -  1  TO  5 
62  Y1  -  F (L) *100  /H 
64  Y  -  100  - Y 1 
66  Z  -  (Y-A) / B 
68  X  -  EXP ( Z) 

70  R(L)  -  (X/B)*(P/H) 

72  NEXT  L 
80  LPRINT  S$, 

82  LPRINT  USING"#####.#  ";  P,H, 

84  LPRINT  USING"####.#  " ; R( 1 ) , R(2 ) ,R(3 ) ,R(4 ) ,R(5 ) 

90  PRINT" ANOTHER  CASE?  Y?  N?" 

92  INPUT  A$ 

94  IF  ( A$  -  "N")  THEN  GOTO  100 
96  GOTO  44 
100  END 


SAMPLE  OUTPUT:  FREQBOUB 


RAIN  RATE(MM/HR)  EQUALLED  OR  EXCEEDED 
Z  HOURS  INDICATED 


STATION 

P 

H 

.1Z 

.  2Z 

.  5Z 

1Z 

2Z 

THULE 

139.7 

211.0 

2.5 

2.1 

1.2 

0.4 

0.1 

PARIS 

569.0 

713.0 

3.5 

3.3 

2.8 

2.1 

1.2 

N. PLATTE 

541.0 

364.0 

6.2 

5.5 

3.9 

2.2 

0.7 

VASH.D.C. 

1064.3 

635.0 

7.4 

6.9 

5.7 

4.1 

2.1 

GUAM 

2260.6 

1177.0 

8.7 

8.4 

7.5 

6.3 

4.4 

CEYLON 

2286.0 

607  .0 

16.5 

15.4 

12.5 

8.9 

4.5 

OIS  RAIN  KATE ( MM/ HE)  EQUALLED  OR  EXCEEDED 
Z  HOURS/TEAR  INDICATED 


STATION 

P 

TRULE 

139 .7 

PARIS 

569.0 

N .PLATTE 

SAI.0 

VASH.D.C. 

1064. 3 

GUAM 

2260.6 

CRTLON 

2266.0 

1 

.1Z 

.  2Z 

211.0 

1.9 

1.3 

713.0 

6.0 

3.9 

364.0 

6.6 

5.1 

635.0 

12.0 

7.2 

1177.0 

14.0 

12.0 

607.0 

28.0 

20.0 

.  5Z 

1Z 

2  X 

0.8 

0.5 

0.3 

2.6 

1.7 

1.1 

3.1 

1.5 

0.7 

5.1 

3.5 

1.7 

8.0 

6.2 

4.4 

13.0 

7.4 

3.7 
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4)  TESTDAY 


2  REM" TESTDAY"  MM/DAY  BY  R.L.WEXLER  NOV  85 
4  OPTION  BASE  I 
6  DATA  6.35,  12.70,  25.4 
8  DATA  3.6,  21.0 

10  DATA  PENSACOLA-WI.297.9,  30.2 

12  DATA  12.8,  8.0,  3.0 

13  DIM  YO ( 5  )  ,  Y ( 5  )  ,  Yl(5) 

14  DIM  F0(5),F(5),F1(5) 

17  FOR  L  -  1  TO  3 

18  READ  R(L) 

20  NEXT  L 

22  READ  A , B 
2A  READ  S$,  P,  D 

25  I  -  P/D 

26  FOR  L  -  1  TO  3 
28  READ  FO (L ) 

30  NEXT  L 

31  LPRINT 

32  LPRINT” TABLE  EXPLICIT  VS  GENERAL  MODEL:  DAILY  RAIN" 

33  LPRINT 

3A  LPRINT  S $ ; " :  P  «"; 

35  LPRINT  USING"#####. #";P; 

36  LPRINT"  MM  D 

37  LPRINT  USING"####. #";D; 

38  LPRINT"  P/D 

39  LPRINT  USING  "###.##";I 
^0  LPRINT 

42  LPRINT"  ","  NO.  DAYS  WITH  RAIN  >  RATE  GIVEN" 

43  LPRINT 

4A  LPRINT" MM/ DAY " , “  OBS","  EST-X" , "  EST-G" 

45  LPRINT 

46  Al  -  38.8-  1.8  *  I 
48  B 1  -  21 

50  FOR  L  -  1  TO  3 

51  Z1  -LOG(R(L)) 

52  Y1  (I. )  -  Al  +  B1*Z1 

53  IF  Y 1 ( L ) <0  THEN  Y1(L)  -0:  GOTO  56 

54  IF  Y 1 ( L  )  > 100  THEN  Yl(L)  -100 

55  FI (L)  -  (10C-Y1 (L))*D/100 

56  NEXT  L 

58  FOR  L  -  1  TO  3 

60  S  -  R(L)/I 

61  X  -  S*P 

62  Y(L )  »  A  +  B  *  (LOG(X)) 

63  IF  Y (L ) <0  THEN  Y(L)  -0:  GOTO  68 

64  IF  Y(L  )>100  THEN  Y(L )  -100 
66  F(L )  -  (100-Y(L))*  D/100 
68  NEXT  L 
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70  FOE  L  -  1  TO  3 
72  LPRINT  USING"####.## 

74  LPRINT  USING"####.# 

76  NEXT  L 

78  LPRINT 

90  FOR  L  -  1  TO  3 

92  YO(L)  "(D-FO (L) ) *IOO/D 

94  NEXT  L 

98  LPRINT 

100  LPRINT  *  "TAB(26)“CUM 

101  LPRINT 

102  FOR  L  -  1  TO  3 

103  LPRINT  USING"####.## 

104  LPRINT  USING"####.# 
108  NEXT  L 

110  LPRINT 

114  LPRINT"  A  B" 

116  LPRINT  USING"##.# 

120  END 


SAMPLE  OUTPUT  X  TESTDAY 


"  ;  FO (L) , FI  (L)  ,  F(L) 


X  RAIN  DAYS" 


-;YO(L),  Y1(L),  Y(L) 


"iAl.Bl 


EXPLICIT  VS  GENERAL  MODEL s  DAILY  RAIN 
a.  PENSACOLA-WI:  P  -  297.9  MM  D  •  30.2  P/D  -  9.14 


MM/DAY 

6.35 

12.70 

25.40 


NO.  DAYS  WITH  RAIN  >  RATE  GIVEN 


OBS 

EST-X 

EST-C 

12.8 

12.1 

12.6 

8.0 

7.7 

8.2 

3.0 

3.3 

3.8 

CUM  2  RAIN  DAYS 


6.35 

12.70 

25.40 


57.6 

73.5 

90.1 


59.9  58.3 
74.4  72.8 
89.0  87.4 


Al  Bl 

21.0  21.0 
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SAMPLE  OUTPUT:  TESTDAY 


EXPLICIT  VS  GENERAL  MODEL:  DAILY  RAIN 

b.SAN  JUAN-YR:  P  -  1499.9  MM  D  -  199.2  P/D  -  7.53 

NO.  DAYS  WITH  RAIN  >  RATE  GIVEN 

MM/DAY  OBS  EST-X  EST-G 

6«35  64.8  71.6  71.8 

12.70  34.0  42.6  42.8 

25*40  15*0  13.6  13.8 


CUM  Z  RAIN  DAYS 

,5*35  67.5  64.1  64.0 

12.70  82.9  78.6  78.5 

25-*°  92.5  93.2  93.1 

Al  B1 

25.2  21.0 


c.  WASH. D . C  .  -SU :  P 


MM/ DAY 

6.35 

12.70 

25.40 


-  316.7 

MM  D  ■ 

NO.  DAYS 

WITH  RAIN 

OBS 

EST-X 

12.8 

12.2 

8.0 

8.0 

3.0 

3.7 

29.2  P/D  -  10.85 

>  RATE  GIVEN 

EST-G 

12.8 

8.5 

4.3 


CUM  Z  RAIN  DAYS 


6-35  56.2  58.1 

12.70  72.6  72.7 
25.40  89.7  87.2 


56.3 
70.8 

85.4 


Al  B1 

19.3  21.0 
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5)TK8TIOUB 

1  REM  *TESTIOUB*  MM/BI  RLV  JAI  30  86 

2  LPBIIT* SAMPLE  OUTPUT:  TltTUOUB* 

3  LPRIBT 

4  IBM  P/I  2.8  -4.3  MM/ II:  A-64.7-6.4(P/l),  1-22 
6  DATA  1.27, 2. 34, 5. 08 ,4. 33,12.70 

10  DATA  PESHAWAR, 373. 4,  108 

12  DATA  38.3,  38.3,  23.0,  18.4,  3.8 

13  DIM  T(5),Tl(3),  T2<3> 

14  DIM  1(5),  A(3) ,8(3) ,  P(S) ,P1(S) ,P2(5) 

17  POB  L  -  l  TO  5 

18  BEAD  1(L) 

20  IBXT  L 

21  HAD  8|,  P,  I 

22  I  -  P/| 

23  A( 1 )  -  64.7-8.8  •  (p/|) 

24  1(1)  -  22 

26  POB  L  •  1  TO  3 
28  BBAD  P(L) 

28  T(L)  -  (I-  P(L))*100/I 
30  MBIT  L 

33  LPBIIt"KZPL1CXT  N0DIL  <X)  V|  C1BB1Al  MODEL  (C):  BOORLT  BAII* 

34  LPBIMT  8$;*:  P  -*j 

35  LPBIIT  OSIBC ”#####. #";P( 

36  LPBIIT*  MM  I  -*{ 

37  LPBIIT  USIIC*#####*;|; 

38  LPBIIT*  P/I 

38  LPBIIT  US1I6  *###.##*;! 

40  LPBIIT 

42  LPBIIT”  * . ' 10 .  IOUBS  WITH  BAII  >  BA I I -BATS  CIVBI* 

43  LPBIIT 

44  LPBIIT* MM/ IB*,  *  088* ,  *  I8T-**,*  I8T-0* 

45  LPBIIT 

46  A(2 )  -  3.6 

47  1(2)  -  21 

48  A(3)«56.4 
48  8(3)  •  .2 

50  POB  l  -  1  TO  3 

32  Y1(L)  •  A(l)  ♦  8(1)  •  (LOC(B(L) ) ) 

53  IP  Y1(L)>70  TIBI  Y1(L)  •  A(3)*B(L)«8(3 ) 

34  IP  Yl(L)>100  TUI  Yl(L)  *100 

35  PI (L)  »(100-Y1(L))*I/100 
56  IBXT  L 

SB  POB  L  *  1  TO  3 
60  8  -  B(L)/I 

62  X  »  S  *  8(2) 

64  Y2(L)  •  A(2 )  ♦  8(2)  •  (LOC(X)) 

63  IP  Y2(L)>100  TIBI  Y2(L)*100 
66  P2(L)  -  (100  -  Y2(L) )  •  1/100 
68  IBXT  L 
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70  FOR  L  -  1  TO  5 

72  LPRINT  USING"####.##  “;R(L), 

74  LPRINT  USING"####.#  " ; F(L) ,Pi (L) ,P2 (L) 

76  NEXT  L 
78  LPRINT 

90  LPRINT  "  ","  " , "CUM  X  FREQ"  * 

92  LPRINT 

93  FOR  L  -  1  TO  5 

94  LPRINT  USING"####.##  ";R(L), 

96  LPRINT  USING"####.#  Y (L) , Y1 (L) , Y2 (L) 

98  NEXT  L 
100  LPRINT 

102  LPRINT  "A  B  a  b" 

104  LPRINT  USING"##.##  ";A(1),B(1),A(3),B(3) 

106  LPRINT"A,B  for  Equation  (4),  a,b  for  Equation  (3)" 
110  END 


SO 


SAMPLE  OUTPUT :  TESTHOUI 


EXPLICIT  MODEL  (X)  VS  GENERAL  MODEL  (6):  HOURLY  RAIH 
a)  TOWNSVILLE:  P  -  1160.8  MM  H  -  364  P/H  -  3.19 

NO.  HOURS  WITH  RAIN  >  RAIN-RATE  GIVEN 


MM/ HR 

OBS 

EST-X 

EST-C 

1.27 

189.0 

186  .0 

188.6 

2.54 

116.0 

130.5 

135.6 

5.08 

66.0 

79.8 

82.6 

6.35 

51.1 

66.9 

65.5 

12.70 

19.0 

22.7 

12.5 

CUM  Z  FREQ 

1.27 

48.1 

48.9 

48.2 

2.54 

68.1 

64.2 

62.8 

5.08 

81.9 

78.1 

77.3 

6.35 

86.0 

81.6 

82.0 

12.70 

94.8 

93.8 

96.6 

A 

B  a 

b 

43.65 

22.00  56.40 

0.20 

A,B  for 

Equation  (4), 

a,b  for  Equation 

(5) 

b)  PESHAWAR:  P  •  373.4  MM  H  -  109  P/H  »  3.43 

NO.  HOURS  WITH  RAIN  >  RAIN-RATE  GIVEN 


MM/  HR 

OBS 

EST-X 

EST-C 

1.27 

59.5 

57.4 

58.1 

2.54 

39.3  . 

40.8 

42.2 

5.08 

25.0 

23.9 

26.4 

6.35 

19.4 

20.0 

21.3 

12.70 

5.8 

6.8 

5.4 

CUM  Z  FREQ 

1.27 

45.4 

47.3 

46.7 

2.54 

63.9 

62.6 

61.3 

5.08 

77.1 

78.1 

75.8 

6.35 

82.2 

81.6 

80.5 

12.70 

94.7 

93.8 

95.1 

A 

B  a 

b 

42.09 

22.00  56.40 

0.20 

A,B  for 

Equation  (4), 

a(b  for  Equation 

(5) 
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6)  INSTANT 


2  REN  "INSTANT"  MM/ MIN  R.  L.  WEXLER  6  DEC  85 

4  OPTION  BASE  l 

5  REM  DATA  FROM  OSWB 

6  REM  TO  EST  THE  DECILES  FOR  MAX-DEPTH  1-HR  STORMS 
8  REM  50.8  MM  (2  INCHES)  AND  127  5  INCHES) 

10  REM  EQUATION  Y  -  A  +  B  LN  X,  X-  S*B,  R  -  S(P/M) 
14  DIM  Y(10),  R(2 ) ,  A(2 ) ,  B(2),  P(2).I(2) 

20  LPRINT 

22  LPRINT  "  RAIN  RATE  DECILES  FOR  MAX  1-HR  STORMS" 
24  LPRINT 

26  LPRINT"  MM/ MIN" 

32  LPRINT"CUM  X  FREQ","  0.85","  2.12" 

34  LPRINT 

40  A(l)  -  -35.1 

42  B(l)  -  28.9 

44  A(2)  -  -69.7 

46  B( 2 )  -  35.8 

48  P(l)  -50.8 

50  P(2 )  -  127 

52  1(1)  -P(l)/60 

54  1(2)  -  P(2  )/60 

58  FOR  L  -  1  TO  9 

60  Y(L)  -  10  *  L 

62  FOR  N  -  1  TO  2 

64  Z  -  (Y(L)-A(N))/B(N) 

66  X  -  EXP(Z) 

68  S  -  X/B(N) 

70  R(N)  -  S  *  I (N ) 

72  NEXT  N 

80  LPRINT  USING"####  ";Y(L); 

82  LPRINT  USING"###.##  ";R(1),  R(2) 

90  NEXT  L 

92  REM  IF  R  -  P/M,  X  -  B 

93  Y1  -  -35.1  +  28.9  *  (LOG(28.9)) 

94  Y2  -  -69.7  +  35.8  *  (LOG(35.8)) 

95  Y1  -  100  -Y1 

96  Y2  -  100  -Y2 

98  LPRINT 

99  LPRINT 

110  LPRINT  "X  FREQ  >  P/M", 

112  LPRINT  USING"###.#  ";Y1,  Y2 

120  END 
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SAMPLE  OOTFOT:  IISTART 


RAIM  BATE  DECILES  VOl  HAS  l-EE  STORKS 

MM/ Mil 


COM  S  FREQ  0.85  2.12 

10  0.14  0.53 
20  0.20  0.72 
SO  0.20  0.96 
40  0.S9  1.27 
SO  0.S4  1.67 
40  0.79  2.21 
70  1.11  2. 93 
•0  1.57  3.87 
90  2.22  5.12 


l  FREQ  >  F/M  37.9  41.6 


7)  ONE-SOUR 

1  REM  "MM/ HIM"  CHANGED  TO  "ONE-HOUR" 

2  REM  "ONE-HOUR*  RLW  86 

3  REM  DATA  FROM  CANAL  ZONE  AND  PANAMA  , 

4  OPTION  BASE  I 

18  DIM  Y(10 )  ,  R(10 ) ,  A( 10 ) , B( 10  )  ,  P(10),  1(10) 

20  LPRINT 

22  LPRINT  “ONE-HOUR  EXCESSIVE  RAIN:  OANAL  ZOMB  AND  PANAMA’ 
24  LPRINT 

28  LPRINT"  "  TAB(  2  8  )  "  RAIN  RATE  f’MM/HR)" 

30  LPRINT  1 

32  LPRINT" CUM  Z  FREQ", 

34  LPRINT"  25.4","  50.8",  "  76.2","  101.6","  127.0" 

36  LPRINT 

38  DATA  15,20,  25,  30,  35,  40,  45,  50,  55,  60 

39  DATA  1,2, 3, 4, 5 

40  FOR  L  °  1  TO  10 

41  READ  Y  (L) 

42  Y(L)aY(L)*(l 00/60) 

44  NEXT  L 

46  FOR  N  -  1  TO  5 

48  READ  P (N) 

49  P (N )  -P(N)*25.4 
54  I (N)  -  P(N) 

56  A(N )  -  -11.5  -  1.15  *  I(N) 

58  B(N)  -  25.3  +  .23  *  I(N) 

60  NEXT  N 

62  FOR  L  -  l  TO  10 

63  FOR  N  -  1  TO  5 

64  Z  -  (Y(L)-A(N))/B(N) 

66  X  -  EXP ( Z) 

68  S  -  X/B(N) 

70  R(N)  «*  S  *  I(N) 

72  NEXT  N 

80  LPRINT  USING"####.# 

82  LPRINT  USING"####.# 

90  NEXT  L 
92  LPRINT 
94  LPRINT 

96  LPRINT 

97  LPRINT  "EQUATION  (1)" 

98  LPRINT  ”  A", 

100  LPRINT  USING'####.# 

102  LPRINT  "  B", 

104  LPRINT  USING"####.# 

110  END 


;  Y(L); 

;R(1),R(2),R(3),R(4),R(5) 


" ; A( 1 ) ,  A(2 ) ,  A( 3 ) ,  A(4),  A(5) 
-;B(1),B(2),B(3),B(4),B(5) 
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ONE-HOUR  EXCESSIVE  RAIN:  CANAL  ZONE  AND  PANAMA 
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